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TERS SARKAC SISTEMI 


Kütle merkezi denge noktasının üzerinde bulunan sistemdir. 


Şekil 1: Ters Sarkaç Modeli 


Normalde sarkaç sabitlendiği noktanın tam altında yer çekimi etkisi ile dengede durmaktadır. Bu noktada 
penduluma etkiyen hiçbir tork olmadığı için sistem dengededir. Eğer ki bu noktadan uzaklaşır ise eski denge 
haline dönebilmek için bir onarma torku oluşur. Rijit bir cisim ile destekli sarkaç tam olarak denge 
noktasının üzerinde iken denge pozisyonundan 180 derece ötededir. Bu sarkaç için dengesi zor bir 
pozisyondur. Bu noktada sarkaç üzerinde hiçbir tork olmamasına rağmen ufacık bir uzaklaşma yerçekimsel 
torka sebebiyet vermekte ve bu noktadan ivmelenerek uzaklaşmasına neden olmaktadır. 
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Şekil 2:Ters sarkacın açısal gösterimi 


Bu ters pozisyonda sarkacı dengede tutmak için, pendulumun düşmeye başladığı andaki açısını ve yönünü- 
kullanan bir geri beslemeli kontröl sistemi kullanılabilir. Ters sarkaç dinamik ve kontrol teorisinde yaygın 
olarak kullanılan bir problemdir. 


Birçok kontrol algoritmasını test amaçlı kullanılabilir. (PID, durum uzay gösterimi , neural ağlar, fuzzy 
kontröl, genetik algoritmalar, vb.) 

Ayrıyeten bu problemde çoklu linkler, üzerinde sarkacı tutan ve hareket edebilen bir araba, bir tahteravalli 
üzerinde kayabilen bir arabanın üzerinde duran sarkaç örnek verilebilir. Ters sarkaç sistemi güdümlü füzeler 
ve roketlerle de bağlantılıdır. Benzer bir problemi anlamak için denge arabaları örnek verilebilir. Parmak 
ucunda dik bir şekilde tutulan bir süpürge sopası, elimizin ufacık bir ileri hareketiyle düşmeye başlar. 
Segway robotları da bu mantıkla işler. 


Segway Balancing a stick Inverted pendulum 


Sekil 3: Denge Robotunun Ters Sarkacla Iliskisi 


Harika Bir Ters Sarkaç Örneği İnsan Vücudu 


Dengelenmiş bir ters sarkaç için verilebilecek en iyi örnek insan olabilir. 


center of mass 


pivot 


Şekil 4: İnsan-Sarkaç modeli 


DA 


Ayakları yere denge noktası (pivot noktası) olarak sabitlenmiş, dik bir şekilde duran insan küçük kas ayarları 
yapmaz ise düşecektir. 

Gözlerden, kaslardan ve eklenmlerden alınan pozisyon bilgisini ve de 

iç kulak içerisindeki yarı dairesel kanallardan ve otolit organlardan alınan yönelim(oriyantasyon) bilgisini 
kullanan denge duyusu bizi dik tutmak için iskeletsel kaslara küçük ayarlamalar yapmaktadır. İşte bizi ayakta 
dik tutan da budur. 

Yürüme, koşma, tek ayak üzerinde durma gibi eylemler bu sisteme eklentiler yapmamızı gerektirir. 

Bölgesel bir rahatsızlık, alkol, madde kullanımı reflekslerimize karışarak dik duramamıza sebep olur. 


Sonuçla; ters sarkaç birçok uygulamada var olan bir problemdir. Ve kontrol sistemlerinin 
tasarlanmasında analizinde de test için kullanabileceğimiz bir modeldir. 


Kullanılacak Model 


Bu çalışmada yatayda hareket edebilen bir arabanın üzerine sabitlenmiş ters sarkaç üzerinde durulacaktır. 


Şekil 5: Araba sarkaç Modeli 


Ek bir yardım olmadan, en ufak bir sarsılmayla dengesi bozulur ve düşer. 


Yatay düzlemde araca bağlı olarak sallanan bir sarkacın hareketi iki eksenlidir. Girişi kuvvet (F) ve 
çıkışlarımız sarkacın açısı (0) , aracın pozisyonu (x) olmaktadır. 


M=aracin kütlesi m=sarkacin kütlesini 
Isarkacın uzunluğunu I=sarkacin atalet momentini 
b= araç ve zemin arasındaki sürtünme katsayısını 

tasvir etmektedirler. 


ae 


Sistemin Matematiksel olarak Modellenmesi 


Şekil 6 : Araç sarkaç üzerindeki kuvvet gösterimleri 


Yatayda hareket eden araç için temek hareket denklemi : 


(1) Mi+bi+N =F 
N sarkaç kaynaklı reaksiyon gücünü temsil etmektedir. 


Sarkaç için yatay kuvvetlerin toplamı denklemi : 
(2) N = më + ml0cos0 — ml? sin 0 


2 numaralı denklem 1 de yazılır ise : 


(3) (M + m)z + bz + ml@cos@ — ml@* sind = F 
Bu bizim temel denklemlerimizden biridir. 


Bir diğer denklem için sarkaca yazılabilir. 
(4) P sin + N cos — mgsin@ = ml + më cos 0 


P ve N terimlerinden kurtulmak için sarkacın merkezine göre moment yazılır ise 
(5) Plsin@ — Nlcos@ = 10 

5 ve 4 numaralı denklemlerin kombinasyonu ile : 

(6) (I + ml*)6 4 mal sin 0 = —ml# cos @ 


Bu denklem sistemimizin ikinci önemli denklemidir. 
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Sistemin Matlab desteği ile Modellenip Analizi 


Sistem denklemlerini Matlab ortamında kodlayarak, uygulanacak girişler için vereceği cevapları analiz 
edebiliriz. Bu işlemi “ode45” komutunu kullanarak Matlab içerisinde ve Simulink içerisinde bloklar 
kullanarak gerçekleştirebiliriz. 


“ode45” Fonksiyonu ile Sistemin Modellenmesi 

Sistemimizde arabanın pozisyonu (x) ve sarkacın açısı (0) çıkışlardır. Bu çıkışlar uygulanan kuvvete ve 
zamana bağlı olarak değişmektedir. Sistem deklemlerinde arabanın hızı ve ivmesi, sarkacın açısal hızı ve 
ivmesi de bulunmaktadır. 


x = Ax+ Bu 


y = Cx + Du 


Sistemin zamana ve girişe göre vereceği cevapları analiz edebilmemiz için durum-uzayı gösteriminden 
yardım alabiliriz. 


Bu gösterimde 


arabanın konumu, x(1) 
arabanın hızı , x(2) 
sarkacın açısı, x(3) 


sarkacın açısal hızı, x(4); 
Olarak tanımlayarak bir x Matrisi oluştururuz. 
Ardından ise bu x matrisinin zamana göre türevi olan x matrisini tanımlarız. 


%inverted pendulum state-space matrix derivated by time 
%Xp=[xdot; xdotdot; theta_dot; theta_dotdot] 

%x=pos of cart 

%theta=pos of pend 

function Xp= invert_pend_last(t,x) 


%System Constant 

mc=0.5; %mass of cart 

mp=0.2; %mass of pendulum 

b=0.1; %friction coefficient between ground and cart 
1=0.3; %lenght of pendulum 

I=0.006; %inertia coefficient 

g=9.81; %gravity 

F=1; %Applied Force 

a=mp*1; K= (mpt+mc)*(I+a*1); 


Xp= zeros(4,1); %State-Space Matrix Derivated by time assigmention 

Xp(1)=x(2); 

Xp(2)= ((F + a*x(4)*sin(x(3)) -b*x(2))*(Ita*1) + (a42)*g*sin(x(3))) 7 (K-(a42)*cos(x(3))) ; 
Xp(3)= x(4); 


Xp(4) = (a*cos(x(3))*(b*x(2) - F) - a*sin(x(3))*(a*x(4)*cos(x(3)) + g*mp + g*mc)) / (K - 
(an2)*(cos(x(3)))42); 


end 


Yukarıda verilen kod X matrisinin fonksiyonudur. Göreleceği üzere Xp ile tanımlanan bu matris içinde 
çok sayıda diferansiyel ifade barındırır. Bu matrisi “ode45” fonksiyonu ile çözdürebiliriz. 


5. 


clear all; close all; clc; 

sode45 fonksiyonu ile zamana göre türevli durum-uzay matrisinin çözümü 
50-10 saniye arasında, x-0.5 ve ?=pi/3 başlangıç değerleri için 

[t,x]= ode45(8invert pend last, [0:0.01:10],[0.5, 0, pi/3, 0]); 


figure (1) 


subplot (2,2,1); 


plot(t,x(:,1)) % araç pozisyonunun çizimi 
title('x') 
grid 


subplot(2,2,2); 


plot(t,x(:,2)) % araç hızının çizimi 
title(”x det) 
grid 


subplot(2,2,3); 


plot(t,x(:,3)) % sarkaç açısının çizimi 
title ("theta") 
grid 


subplot (2,2,4); 


plot(t,x(:,4)) % sarkaç açısal hızının çizimi 
title(“theta det”) 
grid 
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Şekil 7 : ode45 fonksiyonu kullanılarak çözülen sistemin grafikleri 


Grafiklere bakarak aracın kuvvet etkisi ile ivmeli bir şekilde pozisyonunun değiştiği, sarkacın ise salınım 
yaptığını söyleyebiliriz. 
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Simulink Yardımı ile sistemin Modellenip Analizi 


66,29 


Sistemimizin çıkış değerleri olan araç pozisyonu “x” ve sarkaç açısı olan “0” değerini elde etmek için 
Simulinkte “Matlab Function” bloğunu kullanabiliriz. 
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Şekil 8: Sistemin Simulink Blok Modeli 


Bu bloğu bu adı üzerinde bir fonksiyon bloğudur. Ve istediğim çıkış değerlerini vermesi için giriş 
değerlerinden oluşan bir fonksiyon tanımlamamız gerekir. 


function [ddx_o,ddT_o]=fcn(f,T, dT, ddT, dx, ddx) 


5System Constant 
mc=0.5; %Smass of cart 
mp=0.2; %Smass of pendulum 


b=0.1; friction coefficient between ground and cart 
1=0.3; Slenght of pendulum 

1-0.006; Sinertia coefficient 

g=9.81; gravity 


ddx o=(f- mp*l*ddT*cos(T) -b*dx + mp*1* (dT*2)*sin(T))/(mctmp) ; 


ddT_o=(-mp*g*1*sin(T) -mp*1*ddx*cos (T) ) /(I + mp*1%2); 


Bu fonksiyondan x dotdot ve T_dotdot(@_dotdot) ifadelerini elde etmekteyiz. Sonrasında bu 
ifadeleri integratörlere sokarak x_dot,x ve T_dot, T değerlerini elde etmekteyiz. 


ae 


f=1 girişi ve 9=pi/3, x=0.5 başlangıç değerleri için sistemin 10 saniyelik simulasyonu yapılırsa 


Şekil 9: Araç pozisyonu grafiği 


Şekil 10: Sarkacın açı grafiği 


Sistemin denklemleri ve başlangıç değerleri simulink ve matlab ortaminda(ode45 için) birebir aynı 
tasarlanmıştır. Hareketlerin fenotipleri aynı olmasına simulink ve matlab ortamındaki grafiklerde farklı 
değerler alınmıştır. Ode45 fonksiyonunun sistemi Simulinkten farklı çözümlemesi olabilir. 
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Sistemin VRLM araciligiyla modellenip animasyonu 


VR Builder eklentisi kullanarak, oluşturulan sistem modellerinin hareketlerini görselleştirebiliriz. Bu isi 
gerçekleştirmek için sistemin 3 boyutlu modeline ve bu sistemin zamana göre çıkışlarını alacağımız simulink 


modeline ihtiyacımız olmaktadır. 


Şekil 11: Ters-Sarkaç Sisteminin VRLM ortamında modellenmiş hali 
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Sekil 12: VR Sink destekli simulink modeli 
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Vr sink modülüne eklenecek olan; 

*“cart translation” (araç taşınımı) için multiplexer (birleştirici) modül tanımlanmıştır. Araca ait 3 eksen 
bulunmaktadır. Aracın -z ekseninde hareketsiz, -y ekseninde ise platform üzerinde sabit durması için 0.6 
birim, -x ekseninde ise hareket etmesi için sistem çıkışlarından x bağlanmıştır. 

*“pend rotation” (sarkaç dönüşü) parametrelerini ise sistemin Theta çıkışından almaktayız. Multiplexer’a 
bağlı olan [0 0 1] matrisi sarkacın, -z eksenine oturtularak döneceğini belirtir. 

Sistem Simulink ortamında 0 - 10 sn arasında ve 0.01 saniyelik açıklıkta simule edilince, VR animasyonu 


ekil 13: Başlangıç Pozisyonu Şekil 14 : Pozisyon 1 
S y y 
Şekil 17: Pozisyon 2 Şekil 16 : Pozisyon 3 
Şekil 17: Pozisyon 4 Şekil 18 : Pozisyon 5 


aw, 


Başlangıçtan itibaren modelin hareketi incelendiğinde, aracın sağa doğru yöneldiği ve sarkacın salınım 
yaptığı gözlenmektedir. 


GUI 


Matlab GUI sistemlerin incelenmesini kolaylaştırabileceğimiz bir arayüz eklentisidir. 


Şekil 18:Sistemin GUI ortamında modellenmesi 
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Tag: figurel Current Point: [780, 193] Position: [680, 419, 957, 679] 


Modellenen bu sistem sistemin simulink devresiyle etkileşimdedir. Bu etkileşimi GUP nin içine gömülen 
komut ile belirttik. 

Bu arayüzden istediğimiz şudur: 

İlk olarak çıktısını almak istediğimiz değerleri ‘radio buton’lardan seçmekteyiz. Ardından 

Simulation butonuna basıldığı zaman ‘edit textbox’ lar ile belirtilen kutucuklara girilen değerleri ve ‘slider’ 
ile belirtilen simulasyon zamanını simulink sistemine gönderilerek simulasyon gerçekleşmektedir. İstenilen 
grafikler axes2 isimli grafikte yer alacaktır. 
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Şekil 19: Verilen girişler için X ve 0 çıktısı 
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Şekil 20: Verilen 
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SONUÇ 


Ters Sarkaç sisteminin kontrolsüz durumda uygulanacak kuvvete vereceği cevaplar Matlab ortamında analiz 
edilmiştir. Bu analiz sonuçlarına göre sistemin nasıl bir kontrole ihtiyaç duyulacağını kestirebiliriz. 

Mesela sürekli salınımlar yapan sarkacın sönümlenmesi gerekir. Bu sönümleme işlemini aracın konumu 
değiştirerek yapabiliriz. Bunun için sisteme uygulanacak kuvvetin yönü ve büyüklüğü belirlenir. 

Bu sistemden hareketle, başka sistemler içinde kontrol algoritmaları geliştirebiliriz. 
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